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Abstract 

The solution of a (stochastic) differential equation (SDE) can be lo- 
cally approximated by a stochastic expansion, a linear combination of 
iterated integrals. Quantities of interest, like moments, can then be ap- 
proximated with the expansion. We present a formula for the case where 
the drivers of the equation are time and Wiener processes. We also present 
a Mathematica implementation of the result. 

1 Introduction 

We consider the stochastic differential equation 

n-i 

dY t = f{Y t )dt + Y J 9{Y t )dWi 



where the drivers W % are Wiener processes. By convention, we set time as 
the zeroth driver: dt = dW®. We call J a (t) the Stratonovich iterated integral 
according to the multi-index a, on running from to n — 1. We call I a (t) the 
Ito iterated integral according to the multi- index a: 

I a (t) = /.../ dW£...dW% 

J J0<u 1 <...<u k <t 

J a (t) = /.../ °dW£...odW% 

J J<3<ui<...<u k <t 

with \a\ — I. 

Our aim is to arrive at a formula for the expectation of J a for all a. We 
will use two well-known results: the expectation of Ito iterated integrals I a is 
very easy to calculate and a recursive formula link the two iterated integrals. 
Specifically: 
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1. For Wiener processes, EI a — if 3i such that a* ^ 0. Otherwise, EI a (t) = 
t e /£l 

2. J a = j J a -dW a * + \x(pn-i =at^0)J J a --ds (Lemma 3, BQ) 
From this, we can derive a closed formula for the expectation of any J a . 

2 Formula for EJ a 

Lemma 1. J a = ^cplp + p Q 7o,o,...,o where the first sum is done over some 
words /3 which all contain at least 1 non-zero letter, cp and p a are reals, p a can 
be zero. We denote q a the number of zeros in the last ltd integral. 

All quantities depend on a but the exact values of (3 and cp will be of no 
interest for us here. In other words, the lemma states that J a can be written 
as linear combination of ltd iterated integrals, with at most one of them of the 
form /o,....o- The proof comes easily by recursion on the length of a, allowing 
p a to be equal to 0. This result holds for all type of drivers; it is not limited to 
Wiener processes. 

Proof. The lemma is true for £ < 2: 

Jo = h 

Ji = h 

Jo.o = j Iods + = I ,o 

Ji.o = Jhds + = Zx,0 

J ,i = JI o dW 1 +0 = I 0l i 

Ju = JhdW 1 + ^Jlda = A,i + |/ 

Assume the lemma is true for words up to length < £ and consider a word 
ol of length I. Three cases are possible: 

• at = We have: 

Ja = J J a -ds + 

= JC^2 c pJf3 +Pa-h,...fi)ds 
= y^Cffi^p + p a - Jq,...,Q,Q 

where all /3's have at least 1 non-zero letter. 

• ae-i ^ on ^ We have: 

Ja = J J a -dW ai +0 

= /(£ + P a -h,...fi)dW ae 

= y^C^/ff^qj, + Pa-Iq,...,0,a e 
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• ai-i = ag ^ We have: 

J a = I J a -dW ai + I I J„__ds 



C^cplfj +p a -h.,...fl)dW ai + ~ / (E c ^'^' +Pa— ro,...,o)da 



0,...,0,0 



□ 



Lemma 2. For Wiener processes, EJ a (t) = p a .t qol /q a l or 0. 



This follows directly from lemma Q] by taking the expectation on both sides. 
EJ a ^ iff p a ^ 0, i.e. iff J Q has an I component. 

Lemma 3. For Wiener processes, EJ a ^ iff a is a sequence of either or 
pairs mm, where m is a driver 0. 

Proof. For J a to have an /o,...,o component, it is necessary to have: 

a) cxi = and J Q _ to have an Iq....,o component. 

b) ct£ and a^^i = cxg and Jq, to have an Io,....o component. 

Note that a) and b) are mutually exclusive. □ 

Lemma 4. Whenp a ^ 0, p a = 1/2 #{ I* and q a = (#{cti ^ 0})/2 + (#{ai = 
0}) 

Proof. This comes from the proof for result 3: if we are in case a), we integrate 

J a - with respect to time, so q a = q a - + 1. In case b), we integrate J a with 

respect to time and divide by 2, so q a = q a + 1 and p a = p Q /2. □ 

Theorem 1. Given a word a and assuming Wiener processes: 

f if a is not a sequence of and pairs mm 
a \ ' y p a as defined in result 4 otherwise. 

Thus, for example: 

£Jo,i,i,o,o = l/2 2 / 2 <( 3 + 2 / 2 V(3 + 2/2)! - ±£ 

■E'-'o, 1,1,0,0,1 = 

£42,1.1,3.3 = l/2 6 /¥°+ 6 / 2 )/(0 + 6/2)! - ^ 

£^2,2,0,1,1,3,3,0,0,0 = l/2 6 / 2 <( 4 + 6 / 2 V(4 + 6/2)! = ^ 
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3 Mathematica implementation 



The result in theorem Q] can easily be implemented in Mathematica. This iter- 
ative implementation is to replace a recursive one, as it requires less memory 
and is faster, computing the expectation in at most £ iterations (where t is the 
size of word a). 

ExpSBM[t_, j[a_List]] := Module [{i, c}, 
i = LengthOa; 
c = {0, 0}; 
Catch [ 
While [i > 0, 
If [a[[i]] == 0, 
c += {0, 1}; i— , 

IfECi > 1) && (a[[i]] == a[[i - 1]]), 
c += {1, 1}; i -= 2, 
c = {Infinity, 0}; ThrowOO 
]]]]; 

(1/2) ~First@c t~Last@c/(Last@c) !] ; 

Using this code yields the expected results: 

{ExpSBMCt, j[{0, 1, 1, 0, 0}]], 

ExpSBM[t, j[{0, 1, 1, 0, 0, 1}]], 

ExpSBMCt, j[{2, 2, 1, 1, 3, 3}]], 

ExpSBMCt, j[{2, 2, 0, 1, 1, 3, 3, 0, 0, 0}]]} 

{t~4/48, 0, t'3/48, t-7/40320} 

4 Conclusion 

Using a relation between Stratonovich and Ito integrals, we derived a simple 
formula for the expectation of Stratonovich iterated integrals when the drivers 
are Wiener processes. This result was then implemented in Mathematica. 
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